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Abstract 



We propose a new scenario for compound chondrule formation 
named as "fragment-collision model," in the framework of the shock- 
wave heating model. A molten cm-sized dust particle (parent) is dis- 
rupted in the high-velocity gas flow. The extracted fragments (ejec- 
tors) are scattered behind the parent and the mutual collisions be- 
tween them will occur. Wc modeled the disruption event by analytic 
considerations in order to estimate the probability of the mutual col- 
lisions assuming that all ejectors have the same radius. In the typical 
case, the molten thin (~ 1mm) layer of the parent surface will be 
stripped by the gas flow. The stripped layer is divided into about 200 
molten ejectors (assuming that the radius of ejectors is 300 ^m) and 
then they are blown away by the gas flow in a short period of time 
(~ 0.01 s). The stripped layer is leaving from the parent with the ve- 
locity of ~ 4cms~^ depending on the viscosity, and we assumed that 
the extracted ejectors have a random velocity Av of the same order of 
magnitude. Using above values, we can estimate the number density 
of ejectors behind the parent as Ue ~ 800 cm~^. These ejectors occupy 
^ 9% of the space behind the parent in volume. Considering that the 
collision rate (number of collisions per unit time experienced by an 
ejector) is given by .Rcoii = <7coii^^e^^^) where cJcoU is the cross-section 
of collision (e.g., Gooding Sz Keil 1981, Meteoritics 16, 17), we obtain 
-RcoU ~ 36 collisions/s by substituting above values. Since most colli- 
sions occur within the short duration (~ 0.01 s) before the ejectors are 
blown away, we obtain the collision probability of Pcoii ~ 0.36, which 
is the probability of collisions experienced by an ejector in one dis- 
ruption event. The estimated collision probability is about one order 
of magnitude larger than the observed fraction of compound chon- 
drules. In addition, the model predictions are qualitatively consistent 
with other observational data (oxygen isotopic composition, textural 
types, and size ratios of constituents). Based on these results, we con- 
cluded that this new model can be one of the strongest candidates for 
the compound chondrule formation. 

It should be noted that all collisions do not necessarily lead to the 
compound chondrule formation. The formation efficiency and the fu- 
ture works which should be investigated in the forthcoming paper are 
also discussed. 

Keywords: meteorites. Solar System origin. Solar Nebula 
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1 Introduction 



Chondrules are millimeter-sized, once-molten, spherical-shaped grains 
mainly composed of silicate material. They are abundant in chon- 
dritic meteorites, which are the majority of meteorites falling onto the 
Earth. They arc considered to have formed from chondrulc precursor 
dust particles about 4.56 x 10^ yr ago in the solar nebula (Amelin et al. 
2002); they were heated and melted through flash heating events in 
the solar nebula and cooled again to solidify in a short period of time 
(e.g., Jones et al. 2000, and references therein). Typical chondrules 
are single spherical objects, while compound chondrules are composed 
of two or more chondrules fused together. They are rare in all chon- 
drules (~ 4%), but occur in many classes of chondrites, so they offer 
crucial information regarding the physical and chemical state of solid 
materials during chondrule formation (Gooding & Keil 1981, Wasson 
et al. 1995, Sekiya & Nakamura 1996, Ciesla et al. 2004, Akaki & 
Nakamura 2005). 

Wasson ct al. (1995) defined compound chondrules based on tex- 
tures and assigned each constituent chondrule as primary or secondary. 
The primary chondrule was rigid enough to retain its original shape 
at the time of compound chondrule formation. On the contrary, the 
secondary chondrule had a low viscosity enough to allow it to conform 
the shape of the primary. Akaki &; Nakamura (2005) considered four 
types of compound chondrules as follows: (1) the enveloping type, 
wherein a secondary chondrule encloses a primary one, (2) the adher- 
ing type, wherein a small secondary chondrule forms a hemispherical 
bump on the surface of a larger primary chondrule, (3) the consorting 
type, wherein the conjoined chondrules are of similar size, and (4) the 
blurred type, wherein compound chondrules have a blurred boundary 
between two chondrules, which makes it difficult to define the primary 
or the secondary. They classified 30 compound chondrule sets found 
in two CV3 chondrites into 7 blurred, 13 adhering, 4 consorting, and 
6 enveloping. 

At least four models have been proposed for compound chondrule 
formation. (1) Random collision model; totally or partially molten 
particles randomly collided in the solar nebula (Gooding & Keil 1981, 
Sekiya and Nakamura 1996). (2) Collision between individual chon- 
drules in shock- wave heating (Ciesla 2006). (3) Eruption model; im- 
mediately after the formation of the solidified shell in the primary 
chondrule, the inner residual melt was vacuumed out through a crack 
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of the shell, and formed the secondary chondrule on the surface of the 
primary (Sanders &; Hill 1994). (4) Relict grain model; fine-grained 
dust particles accreted on the surface of already-formed primary were 
melted in second heating event (Wasson 1993, Wasson et al. 1995). 
The relict grain model is good for enveloping type if considering that 
fine-grains accreted uniformly on the primary surface, but it seems to 
be difficult to explain other types. In the eruption model, the sec- 
ondary was formed from the inner residual melt after the formation of 
the solidified shell on the primary surface, however, the chemical com- 
positions of some sets of compounds are inconsistent with this scenario 
(Akaki &; Nakamura 2005). The random collision model is considered 
not to account for the observed fraction of compound chondrules be- 
cause of the low density of matter in the nebula (Gooding &; Keil 1981, 
Wasson ct al. 1995, Sekiya and Nakamura 1996). Collision between 
individual chondrules in shock-wave heating can account for the ob- 
served fraction of compound chondrules if the enhancement of the dust 
particle in the pre-shock region is ~ 400 times that expected under 
canonical conditions (Ciesla 2006). However, such a highly dust-rich 
region is unfavorable to explain the scarcity of isotopic fractionation 
of sulfur if the shock wave has a large spatial extent in solar nebula 
(Tachibana & Huss 2005, Miura & Nakamoto 2006). 

In this paper, we propose a new scenario for compound chondrule 
formation. The shock-wave heating model is one of the most plausible 
models for chondrule formation (e.g., Connolly k. Love 1998). In this 
model, the dust particles are exposed to a high-velocity gas fiow and 
heated by the gas frictional heating. It has been suggested that the 
maximum sizes of chondrules are regulated by the gas flow because 
large dust particles should be disrupted by the strong gas ram pres- 
sure when they melt (Susa k, Nakamoto 2002, Kato et al. 2006). Re- 
cently, we carried out three-dimensional hydrodynamics simulations of 
molten dust particle exposed to the gas flow and showed that molten 
cm-sized dust particle is disrupted into many small pieces in a typi- 
cal setting of nebula shocks (Miura k, Nakamoto 2007) . These pieces 
have many chances of mutual collisions to form compound chondrules 
because the local number density of them behind the disrupted par- 
ticle is enhanced. We name this scenario "fragment-collision model" 
and think that it can be a strong candidate for compound chondrule 
formation model. This model seems very similar to the model of col- 
lision between individual chondrules in the shock-wave heating at the 
point that compound chondrules are formed in the gas flow (Ciesla 
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2006). The difference is that in our model, compound chondrules are 
formed from a single large dust particle, so the dust enhancement in 
the pre-shock region is not necessarily required. 

The purpose of this paper is to estimate the collision probability, 
i.e., the number of collisions which a single piece will experience in a 
disruption event, by using a simple formulation. We call the disrupted 
dust particle as "parent" and small pieces as "ejectors." For simplicity, 
we assume that all ejectors have the same radius of in this paper. We 
describe the formulations for estimating the collision probability in ^ 
The expected probability is discussed in ^ We consider appropriate 
situations in which our model can be applied in ^ We compare our 
model with observational data of compound chondrules in ^ Finally, 
we make conclusions in ^ 

2 Formulation 
2.1 Collision Rate 

The number of collisions per unit time experienced by each ejector 
(collision rate) is given by -RcoU = fcoU'T'eAf , where fJcoU is the colli- 
sional cross-section (<Tco1i = 47rrg), ng is the number density, and 
is the velocity dispersion of ejectors (Gooding & Keil 1981, Sekiya & 
Nakamura 1996). The original point of our model is to estimate Ue 
resulting from disruption of the parent. Considering the total number 
of ejectors torn away from the parent in an infinitesimal duration 5t, 
6N, and the volume of the region in which these ejectors are scat- 
tered, 6V, we obtain ric = 5N/5V . Assuming that all ejectors just 
after ejection are parting from the parent with a velocity of ~ Aw, we 
obtain the volume in this phase as 5Vq ~ 2-Kr'^/\v5t^ where rp is the 
radius of parent (see Fig. [1^). It should be noted that ejectors are 
jumping out of rear side of the parent surface, not of front one. After 
ejection, the motions of ejectors are affected by the ambient gas flow. 
We simply assume that ejectors are accelerated with a constant accel- 
eration a in the direction of the gas flow (z-axis), on the other hand, 
in the direction perpendicular to the gas flow (r-axis) they move with 
a constant velocity of ~ At; (see Fig. [Da). The acceleration is given by 
a = 3pfm/4re/5mat5 where is the gas ram pressure and /Omat is the 
material density inside of the molten dust particle. In this later phase, 
the region in which ejectors are scattered is getting wider steeply with 
time t and its volume is given by 5Vt ~ 7r(Ai;t)^at5t. Approximating 
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6V ~ 6Vo + 6Vt, we obtain the number density of ejectors as 

-1 



Hp, ~ 



(i?eject/27rr2At;) [l + {t/Q 



and the colhsion rate 



i^cou ~ (2r^p Reject /rl) [l + {t/uY 



(1) 



(2) 



where -Reject is the total number of ejectors extracted from the parent 
per unit time (ejection rate) defined by -Reject = 5N/6t. The time 

/ a \l/3 

is defined as = ( 8pmat''p'^e/3At>pfm j and it gives the timescale 
within which most collisions will occur (see ^2.2p . 
[Figure 1] 



2.2 Collision Probability 

The probability of collisions that each ejector experiences during the 
time from ti to ^2 is given by integrating the collision rate -Rcoii over 
the time t from ti to t2 as (see Eq. [2]) 

Pcon{ti,t2)= / Rcoiidt. (3) 
Jh 

In order to simplify the integration, we approximate -RcoU as 

_ j 2r^ Reject /rl for t < t,, 

""^"ii - \ 2rlReject/rl {t/uf for t > U. 

Using above approximation, we obtain -Pcoii(^i; ^2) as 

{2rlReyect/rl {t2 - h) when h < t2 < U, 

rlRejeet/rl 3t* - 2ti - {t2/U)~^t2 when ti < < ^2, 
r^Rejecttl/rl (tf - q^) when t,<ti<t2. 

(5) 

In this paper, we count all collisions that expected to occur behind 
the parent. Finally, we obtain the probability of collisions experienced 
by each ejector in one disruption event (collision probability) as 

D D /n ^ 3r^iiejccti* /'''2Pmat?'I -Reject \ ^ /„x 
^'coll = ^'coU (0, 00) = — ^ = I _ ) . (6) 



7 



In addition, we find that -Pcoii(0, t*)/fcoii = 2/3. It suggests that 
most (67% of) collisions occur within i*, therefore, the time can be 
considered as the typical timescale of the mutual collision^ 



2.3 Velocity Dispersion and Ejection Rate 

In order to estimate Av and i?ojcct in Eq- dS]), we consider the hy- 
drodynamics of molten parent particle exposed to gas flow. Before 
melting, since the cm-sized parent is too large to homogenize internal 
temperature due to the thermal conduction, the temperature is higher 
at the parent surface facing to the gas flow than at the center (Yasuda 
& Nakamoto 2005, 2006). It causes to form liquid layer at the parent 
surface. Kato et al. (2006) obtained the internal velocity of the liquid 
layer by analytically solving the steady hydrodynamics equations for a 
core-mantle structure with a linear approximation. According to their 
results, we can approximate a maximum tangential velocity of the liq- 
uid layer as fmax — 0.112p[^h/ fi, where h is the width of the liquid 
layer and n is the viscosity of molten dust particle (see Appendix [A|) . 
Considering that Av is about the same order of magnitude of ^maxi 
we obtain 

Av ~ Umax ~ Pfm/l/lO/i. (7) 

The total volume of liquid layer can be roughly estimated as ~ hr^. 
When the whole part of the liquid layer fragments into ejectors, the 
total number of ejectors is A'e ~ hr^/r^. Since it is considered that 
the fragmentation proceeds in about a fluid crossing time 

icross ~ rp/i;max ~ lO^rp/pfm/l, (8) 

we obtain the ejection rate as 

i?eject ~ Ne/tcToss ~ TpPfm/l^/lO/ir^ . (9) 



2.4 Width of Liquid Layer 

The liquid layer should be thick enough to cause disruption. Kadono 
& Arakawa (2005) carried out aerodynamic experiments in which a 
liquid layer was attached to solid cores, and the breakup of this layer 



^ There is an exact solution of the integration. The exact solution of Pcou(0;Oo) is 
smaller than that with the approximation of Eq. ^ by about 19%. In addition, the ratio 
-Pcou(0,t,)/Pcoii is about 69%, instead of 2/3 = 67%. 
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occurred by means of the interaction with a high-velocity gas flow. 
They discussed that the breakup did not occur at the Weber number 
defined by = Pimh/js^ where 7s is the surface tension, less than 
10 — 20. This result is similar to the finding that the threshold of 
breakup of liquid droplets without solid cores is ~ 10 at O/^ < 0.1, 
where We = Pfm^p/js is the Weber number for a completely-molten 
particle and Oh = /^/y^/Wt^VTs is the Ohnesorge number (e.g., Fig. 1 
of Hsiang and Faeth 1995). Based on their results, we consider that 
ejection will occur at the time when = 10. Therefore, the width 
of liquid layer is given by 



It should be noted that h does not depend on rp. Strictly speaking, we 
should write h as h ^ min [107s/pfrrn ''^p] because h cannot exceed rp. 
However, we do not consider the case in which /i > rp in this paper, 
so we simply use Eq. (fTOj) . 

2.5 Physical Parameters 

The parent particle is assumed to be mainly composed of forsterite. 
The physical parameters adopted in this paper are pmut = 3gcm~^ 
and 7s = 400 dyne cm~^ (Murase & McBirney 1973). For other physi- 
cal parameters, we adopt rp = 1 cm, rg = 300 ^m, = 3x 10^ dyne cm 
and II = 10^ g cm~^ s~^ as a standard set of parameters. We also show 
results for other sets of parameters in §3.11 

3 Estimation of Collision Probability 

Substituting Eqs. dZ]), ([9]), and (fTO|) to Eq. ([6]), we obtain the expres- 
sion of the collision probability as 



Substituting a standard set of parameters (see ^2.5p . we obtain h ~ 
1.3mm, Av ~ 4.0 cms" \ and Reject ~ 2.0 x 10^ s'^ From Eq. ([I]), 
we obtain the number density of ejectors behind the parent as rig ~ 
796cm-3 for t <C t*, where U ~ 0.013s (see pT]) . The total volume 
of these ejectors is given by (4/3)7rrgne ~ 0.09, so they occupy about 



h ~ 107s /pfm- 



(10) 




(11) 
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9% of the space behind the parent in volume. Finally, the collision 
probability is Pcoii ~ 0.67. Surprisingly, above estimation is larger 
than the observed fraction of compound chondrules by one order of 
magnitude or more. If we assume that all collisions lead to compound 
chondrule formation, this result suggests that the fragment-collision 
model can account for the observed fraction when only ~ 10% or less of 
all chondrules formed via the fragmentation events of cm-sized parent 
dust particles (also see §5.ip . Since the probability is close to unity, 
most of ejectors would experience one mutual collision at least one 
time in average. In addition, some of them might experience multiple 
collisions (more than twice) and form multiple compound chondrules. 

3.1 Parameter Dependence 

The expression of collision probability, Eq. (jlip , is complex to under- 
stand how each parameter affects the result. In order to see in details, 
we explicitly write the dependences of physical parameters as follows: 



h 


OC 




(12) 




oc 




(13) 


Across 


oc 




(14) 


-Reject 


oc 


—3 —1 —1 
'^p'^e Pfm^ ' 


(15) 


SVo 


oc 


^p/^~^ 


(16) 


He 


oc 


^p^^c"^Pfm^ 


(17) 


CcoU 


oc 


rl 


(18) 


Av 


oc 




(19) 


-RcoU 


oc 


^P^^c~VfmV"\ 


(20) 




oc 


^2/3 1/3 -1/3 1/3 
Pfm M ) 


(21) 


-fcoU 


oc 


-1/3-2/3-4/3 2/3 


(22) 



Let us see the dependence on the ejector radius re at first. The total 
number of ejectors contained in the liquid layer, A'e, decreases with 
increase of Tq a,s oc . Th.6 ejection timesccile Across 

does not depend 

on Te. Therefore, the ejection rate i?coii ~ -^e/^cross decreases as re 
increases. We consider only the collision at the region close to the 
parent particle because most (67% of) collisions occur at the early 
phase of i < (see ^2.2p . In this case, since the volume SVq does not 
depend on r^, the number density of ejectors Ue oc Reject /SVo decreases 
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with increase of as oc r^T • The colhsional cross section fXcoU is 
proportional to and the velocity dispersion Av does not depend on 
Tg, so the collision rate i?coii ~ CTcoiin-cAv is inversely proportional to 
Te. We can roughly estimate the collision probability by PcoU ~ ^coU^*) 
which is smaller than Eq. ([6]) by a factor of 2/3, and is proportional 

1/3 

to Te , SO it is found that PcoU decreases with increase of re- 

We will show other dependences of PcoU on the parent radius r^, 
the gas ram pressure pfm, and viscosity of molten parent particle /i as 
below. 

3.1.1 Parent Radius 

Fig. [2] shows the value of PcoU as a function of for various values 
of Tp = 0.2 (dashed), 1.0 (solid), and 5.0 cm (dotted-dashed), respec- 
tively. It is found that PcoU is larger for smaller parent than for larger 
one as estimated by Eq. (j22p . The disruption timescale is longer for 
larger rp, on the other hand, N,. increases with increase of rp more 
steeply, so i?eject increases as rp increases. However, for larger rp, ejec- 
tors are disrupted into more wider region. As a result, rie and .RcoU 
are inversely proportional to rp. The time increases with increase of 
rp, but it cannot completely cancel the dependence of -RcoU- Finally, 
we obtain PcoU oc rp as seen in Eq. ([22l) . For the parameter range 
we adopted in Fig. [2l PcoU ranges about from 0.2 to 2. 
[Figure 2] 

3.1.2 Gas Ram Pressure 

Fig. [3] shows the value of PcoU as a function of for various values 
of Pfm = 10^ (dashed), 3 x lO'' (solid), and 10^ dyne cm"^ (dotted- 
dashed), respectively. It is found that PcoU is larger for weaker gas ram 
pressure than for stronger one as estimated by Eq. (1221) . The stronger 
gas flow can disrupt thiner liquid layer (Eq. [T2|) . so A'e and -Reject 
decreases as pfm increases. This results in that Ue and i?coii are also 
inversely proportional to pf^^- In addition, the strong gas flow scatters 
disrupted ejectors rapidly, so decreases as pfm increases. Finally, we 
obtain PcoU Pim^ as seen in Eq. ([22]) . For the parameter range we 
adopted in Fig. [3l PcoU ranges about from 0.06 to 6. 
[Figure 3] 
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3.1.3 Viscosity 



Fig. m shows the value of i-*coii as a function of for various values 
of /i = 10 (dashed), 10^ (solid), and lO'^ gcm~^ s~^ (dotted-dashed), 
respectively. It is found that PcoU increases as fi decreases as estimated 
by Eq. (122]) . In a highly viscous case, -Reject becomes small because 
of a long ejection timescale (Eq. [T5]) . However, in this case, since the 
velocity dispersion Av also decreases and ejectors are disrupted into 
narrow region, Uq does not depend on fi. Therefore, iZcoU is inversely 
proportional to /i, which is due to the dependence of Av. Considering 
the dependence of t^,, we obtain PcoU oc /i^^/^ as seen in Eq. (p2]) . For 
the parameter range we adopted in Fig. HJ PcoU ranges about from 
0.06 to 6. 

[Figure 4] 



4 Parameter Range 



4.1 Parent Radius 



In Eq. (jlOp . we estimated the layer width h from the disruption con- 
dition of partially-molten parent particle. However, it should be noted 
that this estimation is not valid for much larger parent particle be- 
cause the disruption timescale becomes larger than the thermal con- 
duction timescale. In this case, the layer width h will increase by 
the rapid thermal conduction before disruption as explained below. 
The thermal conduction timescale in the liquid layer tcond is given by 
icond ~ Pm3.tCh?/K ~ 100/9matC7s^/'^Pfm> i^^ sccoud sign of similar- 
ity we substitute Eq. (flO]) . where C is the specific heat per unit mass 
and K is the heat conductivity. On the other hand, we obtain the dis- 
ruption timescale icross ~ ^pfi/^s by substituting Eq. ([TO]) to Eq. ([8]). 
The layer width h would not change significantly during disruption if 
the condition of tcross < icond is satisfied. We rewrite this condition as 

Adopting C = 10^ergg~^K~^ and k = 4 x 10^ ergcm"^ s~^ 
(Murase & McBirney 1973), we find rp ^ 5.3cm for a standard set of 
parameters (see ^2.5|) . If rp is much larger than this value, the layer 
width h changes during disruption because the thermal conduction 
can rapidly transfer energy inside the parent and melt deeply. In this 
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study, we do not treat such situation because we consider cm-sized 
parent particle, which satisfies the condition of Eq. (j23p . 

4.2 Gas Ram Pressure 

The gas ram pressure pfm should be sufficient to disrupt the molten 
parent particle. For simplicity, we consider a completely-molten par- 
ent particle here {h = rp). The disruption will occur when the Weber 
number We = Pfm''p/7s exceeds about 10 (Hsiang & Faeth 1995). 
Therefore, the condition that should satisfy is given by 

Pfm ^ 107s/rp. (24) 
We find ^ 4 X 10'^ dyne cm^^ for a standard set of parameters (see 

Fig. Oshows the gas ram pressure expected to affect the molten 
parent particle just behind the shock front (solid lines) as a function of 
the shock velocity Vs and the pre-shock gas number density uq. The 
gray region is the chondrule-forming shock condition in which the 
gas frictional heating is sufficient to melt the precursor dust particle 
but not so strong as it evaporates the dust completely (lida et al. 
2001). The gas ram pressure is pfm ~ 10^ — 10^ dynecm"^ when no ~ 
10^^ — 10^^ cm~^ and Vg ^ 7 — lOkms^^. The shocks associated with 
gravitational instability (e.g.. Boss & Durisen 2005) or planetesimal 
bow shocks (e.g., Hood 1998, Weidenschilling et al. 1998) might be 
preferable for such situation. 
[Figure 5] 

4.3 Viscosity 

The viscosity of molten parent particle should be small enough to 
cause disruption. The timescale of disruption is given by t 

cross ^ 

fpfJ'hs (see Eqs. [8] and [TO]) . On the other hand, the heating event 
ceases within the timescale that the parent dust particle stops against 
the ambient gas (lida et al. 2001). The timescale is the stopping 
time, which is given by tstop ~ /Omat^'p'yg/Pfmi where is the relative 
velocity between ambient gas and dust particleU If tcross > ^stopj the 
gas drag heating ceases before the disruption occurs. Therefore, the 

^Substituting pfm = Pgi'g, where pg is the gas density, we obtain a well-known expression 
of the stopping time igtop ^ Pmatrp/ PgVg. 
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viscosity at the disruption satisfies the condition of tcross < ^stop- We 
rewrite this condition as 

PmatTs^g ^ (25) 
Pirn 

Substituting Vg = = 8kms~^ (see Fig. [5]), we find fi ^ 3.2 x 
10^ gcm~^ s~^ for a standard set of parameters (see ^2.5p . 

In addition, the viscosity should be large enough to adopt the 
steady solutions (Kato et al. 2006) as discussed below. The liquid 
layer will be ejected within the timescale of icross ~ i~pl^/ls (see Eqs. 
[51 and [TOl) ■ and it decreases as /i decreases. For example, we obtain 
Across ~ 2.5x 10~^ s for Tp = 1 cm and fj, = 0.1 g cm^^ s^^. However, the 
liquid layer cannot be accelerated to the velocity given by Eq. ([7]) in 
such a short duration. Considering that the momentum of the gas flow 
converts to the motion of the liquid layer effectively, the timescale in 
which the liquid layer is accelerated to Av is given by tacc ~ mAv/F, 
where the mass of the liquid layer is m ~ 2TTpaia_trph, the force due to 
the gas flow is F = vrrpPfm, and the velocity dispersion is Av ~ 7s/m 
(see Eqs. 171 and [TO]l . Then we obtain 

tacc ~ 20/9mat7s^/?'fm/^- (26) 

When tacc > Across, the liquid layer does not reach the steady solutions 
before ejection. Therefore, the condition of tacc < Across should be 
satisfied in our model. This condition can be rewritten as 

(27) 

We find fi ^ 2gcm^^ s^^ for a standard set of parameters (see ^2.5|) . 



5 Comparison with Observations 
5.1 Fraction of Compound Chondrules 

We found that the collision probability estimated by our model, -Pcoih 
is close to unity, which is larger than the observed fraction of com- 
pound chondrules by one order of magnitude or more (see If we 
assume that all collisions lead to compound chondrule formation, this 
result suggests that the fragment-collision model can account for the 
observed fraction when about 10% or less of all chondrules formed 
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via the fragmentation events. The fraction of chondrules which have 
undergone the fragmentation event depends on a size distribution of 
precursor dust particles in chondrule-forming region. Nomura et al. 
(2007) solved coagulation equations for various sizes of settling dust 
particles in the solar nebula. They obtained the dust size distributions 
for various positions in the nebula and ages of the nebula. According 
to their results, for example, the dust particles with the radii from 
1mm to 1cm have the size distribution similar to dn oc r~^d{logr), 
where r is the dust radius, in regions near midplane at the distance 
of 1 AU or 10 AU and at the age of 10^ yr (it corresponds to the line 
of the slope —1 in Fig. 4 of their paper). This suggests that there are 
10'^ of mm-sized dust particles per a cm-sized one. Assuming that the 
cm-sized particle is disrupted into mm-sized ejectors in a shock- wave 
heatingH total number of ejectors is estimated as ~ 10^. If we also 
assume that all of them become compound chondrules with two con- 
stituents, about 500 sets of compounds are expected to form. In this 
case, the fraction of compound chondrules is ~ 5% (= 500/10^). This 
is very close to the observed fraction of compound chondrules. 

However, all collisions do not necessarily lead to compound chon- 
drule formation if we consider the situation of mutual collision in more 
detail. We discuss the efficiency of compound chondrule formation in 

5.2 Oxygen Isotopic Composition 

In the fragment-collision model we proposed in this paper, the con- 
stituent chondrules of compounds are likely to have similar compo- 
sitions because ejectors originate from the same parent. Akaki & 
Nakamura (2005) measured the oxygen isotopic compositions for 3 
sets of blurred-type compounds, 6 sets of adhering- or consorting-type 

■^In the shock-wave heating model, all of the equations governing the evolution of the 
precursor dust particle are scaled with the initial precursor radius ro as long as the post- 
shock gas properties (temperature, density, and so forth) are spatially uniform (see §4.3 in 
Miura et al. 2002, Icarus 160, 258). It means that the peak temperature of the precursor 
dust particle does not depend on the radius. The reason of the scaling law is that all of 
the heating rates (due to the gas friction, radiation, and so forth) and cooling rates (due 
to the radiative cooling, latent heat of evaporation, and so forth) are given as the rates 
per unit area and the peak temperature is determined by the balance between heating and 
cooling, so the dependence on the particle radius does not appear. Namely, if mm-sized 
dust particles melt (and chondrules are formed), cm-sized dust particles also melt in the 
same shock condition. 
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compounds, and 2 sets of enveloping- type compounds. It was found 
that in a three- isotopic diagram, all sets of blurred-, adhering-, and 
consorting-types fall in the typical range obtained for single chon- 
drules from the same CVS chondrites. These results suggest that the 
two constituent chondrules of these compounds originated from the 
same dust reservoirs as those single chondrules. These observations 
are consistent with our new model for compound chondrule formation. 

In contrast, in one set of the enveloping-types, the oxygen isotopic 
compositions differ between two constituent chondrules. This result 
might suggest that this enveloping-type compound has not formed 
by our model but by the relict grain model, in which fine-grained 
dust particles accreted on the surface of already-formed primary were 
melted in a second heating event (Wasson 1993, Wasson et al. 1995). 

5.3 Textural Types 

There is a dependence of observed fraction of compound chondrules 
on textural types of component chondrules. Gooding & Keil (1981) 
found that compound chondrules with non-porphyritic pairs are more 
frequent than that with porphyritic pairs. From their thin section 
results, they estimated that 13% of non-porphyritic and 2% of por- 
phyritic chondrules are compound or cratered, which are interpreted 
as products of collisions between plastic chondrules. Akaki & Naka- 
mura (2005) reported the same tendency. It is considered that por- 
phyritic and non-porphyritic textures have been formed from partially- 
and completely-molten dust particles, respectively (Lofgren & Rus- 
sell 1986). It is also naturally considered that a partially- molten dust 
particle is highly- viscous, in contrast, completely- molten one has lower 
viscosity. Based on above assumption, the dependence of the observed 
fraction of compound chondrules on the textural types might reflect 
the dependence of the collision probability on the viscosity of molten 
parent dust particle. Actually, our model predicts that lower viscosity 
results into larger collision probability as shown in Fig. HI This result 
is consistent with the observations. 

In addition, various types of textures are seen in compound chon- 
drules, e.g., P-P, P-NP, and NP-NP pairs, where P and NP stand for 
porphyritic and non-porphyritic textures, respectively (Wasson et al. 
1995). In contrast, our model seems not to account for compound 
chondrules with porphyritic textures because the following proper- 
ties are implied: (a) completely-molten ejectors are extracted from 
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completely-molten parent particle, and (b) temperatures of all ejectors 
are the same at the extraction. However, regarding the point (a), it 
can be considered that the partially-molten parent particle, which in- 
cludes tiny unmelted cores inside, behaves as a fluid if the molten part 
occupies the most volume of the parent. In this case, the partially- 
molten ejectors could be extracted from the partially-molten parent 
and the porphyritic textures might be formed. Regarding the point 
(b), we can consider the case that the temperatures of ejectors are not 
uniform at the extraction. If the initial temperatures of ejectors are 
different, various types of textures can be formed in compound chon- 
drules. In order to verify above hypotheses, we must quantitatively 
investigate the thermal evolutions of the parent dust particle and each 
ejector at disruption. It requires to carry out the three-dimensional 
(thermo-)hydrodynamics simulation and/or aerodynamical disruption 
experiment. 

In addition, Connolly & Hewins (1995) reported that porphyritic 
textures can be reproduced from totally-molten droplets by dust seed- 
ing. Based on their results, even if ejectors extracted from the parent 
are totally-molten, they can obtain porphyritic textures inside them 
with the help of the dust seeding. Considering that such seeded ejec- 
tors collide with others before solidification, compound chondrules in- 
cluding porphyritic textures might be formed. Therefore, the effect of 
the dust seeding in compound chondrule formation is also an impor- 
tant issue that should be investigated in the future. 

5.4 Efficiency of Compound Chondrule For- 
mation 

In ^5.11 we implicitly assumed that all collisions between ejectors make 
compound chondrules, however, it is an upper limit of the probability 
of compound chondrule formation: for example, a collision between 
two totally-molten and high-temperature (low-viscous) ejectors does 
not lead to a compound chondrule because they will fuse into a single 
droplet. In addition, a collision after significant cooling has a difficulty 
in compound chondrule formation because both ejectors will solidify 
completely. Therefore, we have to clarify the collision probability 
during ejectors are in the moderate cooling phase. 

Assuming that ejectors cool due to the radiative cooling, the cool- 
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ing rate is given by 

^cool - — , (28) 

where eemit and Cabs are the emission and absorption coefficient, re- 
spectively, dsB is the Stefan-Boltzmann constant, Tg is the tempera- 
ture of the ejectors, Tj-ad is the effective temperature of the ambient 
radiation field (Desch & Connolly 2002, Miura & Nakamoto 2006), 
and C is the specific heat. The cooling timescale in which the tem- 
perature decreases by AT is given by ~ AT/i?cooi- Here, we assume 
that AT = 100 K and 300 K are required for moderate cooling and 
complete solidification, respectively. The effective radiation tempera- 
ture Trad should be lower than 1273 K in order to prevent the isotopic 
fractionation of sulfur in chondrules (Miura & Nakamoto 2006) . The 
cooling rate is reduced by the term of T^ad; however, it is a minor 
effect in the cooling phase (if we assume T^ad = 1000 K, we obtain 
Tjj'ad/Tg' = 0.15), so we neglect this term for simplicity. Assuming 
that eemit = 1 and Tg = 1600 K, compound chondrules can be formed 
only in the phase from t^ooX = 0.08 s to tsoUd = 0.24 s (moderate cool- 
ing phase). Since the timescale of the mutual collision (t* ~ 0.01 s, see 
§2.2|) is shorter than tcool and tsoiidi we obtain the collision probability 
between tcooi and tgoiid as (see Eq. [5]) 

PU = ^'coll (tcool, t.olid) ^ (^] ' , (29) 

V ''COol / 

where we neglect the term of tsoUd because tgojid ^ ^cooV Comparing 
with Eq. ([6]), we find that P'^^w is smaller than Pcoii by a factor of 
Q ^ P'cow/PcoW = (**/icooi)^/3 ~ 5 X 10"^. The physical meaning 
of Q is the efficiency of compound chondrule formation per a mutual 
collision. It means that almost all (99.5% of) collisions do not lead 
to compound chondrule formation. Taking into account the low effi- 
ciency, the probability of compound chondrule formation is about one 
order of magnitude smaller than the observed fraction of compound 
chondrules, although this estimation is larger than that estimated from 
the random collision model (Gooding & Keil 1981, Sekiya & Nakamura 
1996). 

However, we can consider other possibilities which enhance the 
collision probability. One possibility is that the gas flow is blocked 
by the parent particle and ejectors. In this case, the shadowed re- 
gion in which there is no effect of the gas flow appears behind them. 
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Since the ejectors in the shadowed region are not accelerated by the 
gas flow, the concentration of ejectors will be higher than the case 
neglecting this effect. The second possibility is to consider the colli- 
sions between different-sized ejectors. The relative velocity between 
the different-sized ejectors is given by ~ Av + {v^eil, where |t>rei| is 
the bulk relative velocity between different-sized ejectors due to the 
difference of the acceleration (see ^2.ip . This large relative velocity be- 
tween different-sized ejectors is expected to result into large collision 
probability as long as the collisional disruption does not occur (also 
see §5.5p . The third possibility is to consider that the disruption from 
partially-molten parent particle is possible. In this case, the moderate 
cooling might not be needed, so tcooi can be shorter and it leads to the 
increase of P^^^ (see Eq. [29|) . These possibilities are very important 
to investigate in detail, however, they are beyond the scope of this 
paper. We are planning to investigate these issues and the results will 
appear in the forthcoming papers. 

5.5 Size Ratio of Primary and Secondary 

Wasson et al. (1995) measured median diameters of primaries and 
secondaries in compound chondrules and the ratio of the "diameter" 
of the secondary divided by that of the primary. They found that the 
mean and median of the ratio are about 0.3 and 0.25, respectively. 
These results suggest that compound chondrules with different-sized 
pairs appear more frequently than that with same-sized pairs. On the 
other hand, in our model, if there is a size difference between colliding 
two ejectors moving in the gas flow, it leads to a large relative velocity 
between them because the acceleration a depends on the ejector radius 
(see §2.ip . The large relative velocity will promote more frequent 
collisions as commented in ^5.4[ In contrast, the undesirably large 
relative velocity causes disruption of melt droplets upon impact, which 
is an opposite scenario against compound chondrule formation. Here, 
we discuss the appropriate range of diameter ratio of the two ejectors 
that does not exceed the limit of the critical destruction velocity. 

Ejectors are accelerated by the gas flow with the constant acceler- 
ation a = 3pfni/4rePmat (see §2.1). The velocity relative to the parent 
particle in the z-direction, Vz, and the position from the parent, z, are 
given by 

Vz = at = , (30) 

4rePmat 
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(31) 



SrePmat ' 

where t is the time after extraction. It is considered that the compound 
chondrule formation occurs after the ejector cools moderately (see 
^5.4p . Therefore, we consider the region of 



■2- ^cool — 



cool 



Next, ehminating t from Eqs. (jSOp and (jSip . we obtain 



(32) 



(33) 



Since the velocity Vz depends on the ejector radius re, there is the 
relative velocity v^^x between the large ejector with radius r\ and the 
small ejector with radius at the same position. It is given by 



\2np 



mat 



1/2 



1/2 



1 



(34) 



It is found that Urei increases as z increases. The compound chondrule 
formation can occur only if the condition of Urei < t'dest is satisfied, 
where Vdest is the critical destruction velocity, otherwise ejectors will 
be disrupted upon impact. This condition is rewritten as 



^ < ^dcst 



2nPmat^Lt 
3pfm 



1/2 



(35) 



From Eqs. (j32p and ()35p . the region in which compound chondrules 
can be formed is Zcooi < z < 2:dest- In other words, the necessary 
condition for compound chondrule formation is given by Zcooi < -^desti 
which is rewritten asj 



T\ ^ ( 4n/0mat^dcst ^ ^ 
n V 3pfmtcool 



(36) 



Substituting r\ = 300 //m, = 3 x 10^ dyne cm^^, and t^ooi = 0.08s 
to Eq. (j36p , we obtain the appropriate size ratio of the small chondrule 
to large one as r^/r\ > 0.99 for v^est = lO^cms"^, Vg/ri > 0.91 for 
Vdest = lO^cms"^, Ts/n > 0.44 for v^est = lO^cms"^, and rs/n > 



*We substitute ~ n in Eq. ([32 
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0.28 for Vfiest = lO^cms"^, respectively (see Fig. [6]). Ciesla (2006) 
mentioned that as chondrules cool they can survive collisions with 
one another at velocities up to lO^cms"^ due to viscous dissipation 
in the melt. If so, compound chondrule with the size ratio of about 
^ 0.5 can be formed. However, it does not account for the mean value 
of Ts/n from observations (~ 0.3, Wasson et al. 1995). In order to 
explain the observation, it is required that Wdest ~ lkms~^ or more, 
however, it seems unrealistic. 

In order to overcome this difficulty, we can consider the same possi- 
bilities as discussed in §5.41 For example, in the region shadowed from 
the ambient gas flow, Vrei does not increase with z because ejectors 
are not accelerated. If v^ei just after entering the shadowed region is 
less than Vdest, we obtain z^est oo because v^ei does not increase fur- 
ther. In this case, the compound chondrule formation is possible for 
arbitrary rg/n. It implies that the compound chondrules with smaller 
value of Tg/ri can be formed in the shadowed region. In addition, 
we discussed the possibility that the disruption from partially-molten 
parent. In this case, tcooi can be shorter than that we assumed in this 
subsection. The shorter tcooi allows the compound chondrule forma- 
tion for the wider range of the size ratio rg/ri (see Fig. [6]). These 
issues will be discussed in detail in the forthcoming paper. 
[Figure [6] 

6 Summary 

We proposed a new scenario for compound chondrule formation named 
as "fragment-collision model," in the framework of the shock-wave 
heating model. We modeled the disruption of molten cm-sized par- 
ent dust particle exposed to a high-velocity gas flow in order to esti- 
mate the efficiency of mutual collisions between small fragments as- 
suming that all of them have the same radius. We obtained colli- 
sion probability PcoU for a wide range of parameters (parent radius 
Tp = 0.2 — 5 cm, ejector radius = 100 — 1000 ^m, ram pressure of 
the gas flow = 10^ — 10^ dynecm"^, and viscosity of molten par- 
ent = W — 10^ gcm~^ s^^). The estimated collision probability was 
~ 0.1 — 1 or more, which is about one order of magnitude larger than 
the observed fraction of compound chondrules. If we assume that all 
collisions lead to compound chondrule formation, this result suggests 
that our model can account for the observed fraction of compound 
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chondrules when about 10% or less of all chondmlcs formed via the 
fragmentation events. Since the fraction of chondrules which have 
undergone the fragmentation event depends on a size distribution of 
precursor dust particles in chondrule-forming region, it would be diffi- 
cult to make a conclusion about the fraction of compound chondrules. 
However, numerical results of dust coagulation equations in the solar 
nebula seem to match well with our estimation (Nomura et al. 2007). 
In addition, our model does not require the dust enhancement in the 
pre-shock region because compound chondrules are formed from a sin- 
gle large dust particle. This is advantageous to explain the scarcity 
of isotopic fractionation of sulfur (Tachibana & Huss 2005, Miura & 
Nakamoto 2006). 

We also compared our model with other observational data. Akaki 
& Nakamura (2005) measured the oxygen isotopic compositions of con- 
stituents of compound chondrules and found that in a three-isotopic 
diagram, all compound chondrules except for enveloping-types fall in 
the typical range obtained for single chondrules from the same chon- 
drites. These observations are consistent with our model because two 
constituent chondrules are expected to originate from the same dust 
reservoirs. Gooding & Keil (1981) and Akaki k, Nakamura (2005) re- 
ported that compound chondrules with non-porphyritic pairs are more 
frequent than that with porphyritic pairs. It is thought that these ob- 
servations can be also explained by our model because the collision 
probability depends on the viscosity of molten parent dust particle. 
The dependence on the viscosity is consistent with the experimen- 
tal results that porphyritic and non-porphyritic textures have been 
formed from partially- and completely-molten dust particles, respec- 
tively (Lofgren k, Russell 1986). Finally, the size ratios of secondary 
to primary in each set of compound chondrule have the mean value 
about 0.3 (Wasson et al. 1995). This result might be explained by our 
model because two fragments with different sizes are accelerated by 
the gas flow with different accelerations. As a result, these two frag- 
ments obtain large relative velocity and it would enhance the collision 
probability. Therefore, the compound chondrule of different-sized pair 
tends to be formed more frequently than that of same-sized one. 

However, it should be noted that all collisions do not necessar- 
ily lead to compound chondrule formation. For example, undesirably 
fast collisions cause disruption of ejectors upon impact. Assuming that 
the upper limit of the collisional velocity for coalescence is lO^cms^^ 
(Ciesla 2006), the appropriate size ratio of secondary (small) chon- 
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drule to primary (large) one is rs/n ^ 0.5, which does not account 
for the observations (mean value of rg/n ~ 0.3, Wasson et al. 1995). 
In addition, ejectors should cool moderately before collision to make 
compound chondrule not to fuse into a single droplet. The collision 
probability after the moderate cooling, however, is much smaller than 
that of total collisions. In order to overcome these difficulties, we con- 
sider other physics that we did not take into account in this paper 
(e.g., the gas flow is blocked by the parent and/or numerous num- 
bers of ejectors). These issues will be discussed in the forthcoming 
paper. In addition, our model should be tested by other methods, 
e.g., three-dimensional (thermo-)hydrodynamics simulation or aero- 
dynamic disruption experiment in the future. 
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A Hydrodynamics in Liquid Layer 

Kato et al. (2006) examined the hydrodynamics of the liquid layer by 
analytically solving the hydrodynamics equations for a core-mantle 
structure with a linear approximation. Fig. [7] shows a schematic pic- 
ture of the set-up in their analysis. According to their solutions, we 
can obtain the radial and tangential velocities, Vr and vg, at arbitrary 
position in the liquid layer. The tangential velocity has the maximum 
value Vmax at = 0.2657r on the surface of liquid layer. Fig. [8] shows 
Vmax. as a function of the layer width h. The horizontal axis is a nor- 
malized layer width h/rp and the vertical one is a normalized velocity 
Vmax/ (Pfm^'p/ m)- The solution obtained by Kato et al. (2006) is a com- 
plex function (dashed), however, it can be approximated by a linear 
interpolation between Umax = for h/r^ = and Vmax = 0.112pfjnrp/|U 
for h/vp = 1. The value of fmax for /i/rp = 1 corresponds to the so- 
lution obtained by Sekiya et al. (2003), in which they analyzed the 
hydrodynamics of a completely-molten dust particle. In our model. 
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we adopt the linear interpolation for fmax (solid line) as the velocity 
dispersion of ejectors disrupted from a molten parent particle, which 
is given by 

Vma.^ = 0.112pfrah/ll. (37) 

[Figure \7\ 
[Figure [H] 
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(a) just after ejection (b) far from parent 



gas flow 




Figure 1: Schematic picture of fragmentation of molten chondrule precur- 
sor dust particle (parent) in the shock-wave heating event. When the gas 
ram pressure is too strong for the molten particle to keep its shape, many 
small particles (ejectors) are dispersed behind the parent particle and mu- 
tual collisions between them will occur, (a) Just after ejection, these ejectors 
part from the surface of the parent at rear side with initial ejection veloc- 
ity, (b) Far from the parent particle, the motions of ejectors become an 
uniformly-accelerated motion along the direction of the gas flow (z-axis), 
and a constant- velocity motion perpendicular to the gas flow (r-axis). 
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Figure 2: Collision probability Pcou are plotted as a function of ejector radius 
Te for various parent radii of rp = 0.2, 1.0, and 5.0 cm. We adopt standard 
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values for other parameters, pfn 
(see g23D. 
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Figure 3: Same as Fig. [2] except for various gas ram pressures of pfn 
3 X 10^, and lO^dynecm"^. Other parameters are described in §2.51 
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Figure 5: Gas ram pressure affecting the mohen parent particle behind the 
shoclc front. The horizontal axis is the shock velocity Vg and the vertical axis 
is the pre-shock gas number density uq. The gray region is the chondrule- 
forming shock condition in which the gas frictional heating is sufficient to 
melt the precursor dust particle but not so strong as it evaporates the dust 
completely (lida et al. 2001). 
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Figure 6: Appropriate ranges of size ratio of the small chondrule (radius of 
Ts) to large one (radius of ri) for compound chondrule formation are displayed 
as a function of assumed destruction velocity f dest in the case of ri = 300 /im 
and pfm = 3 X 10"^ dyne cm~^. According to Ciesla (2006), fdest ~ lO^cms"^. 
The solid, dashed, and dotted curves are criteria of the destructive collision 
for the timescales of the moderate cooling tcooi of 0.008 sec, 0.08 sec, and 
0.8 sec, respectively, below which the relative velocity between the small and 
large chondrules is too large for coalescence. The gray region suggests the 
mean or median values from the observational data (Wasson et al. 1995). 



32 



gas flow 



solid 




liquid layer 



Figure 7: Schematic picture of the set-up in Kato et ah (2006) to solve 
the hydrodynamics equations of molten dust particle. A solid core (gray) of 
radius Tc is surrounded by a hquid layer (white) of width where the dust 
radius r-p = rc + h. The radial and tangential fluid velocities inside the liquid 
layer, Vr and vq, were obtained as a function of r and 6, where r is a distance 
from the center. The ambient gas flows from up to bottom in this figure. 
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Figure 8: Maximum tangential velocity at the surface of liquid layer Vmax as 
a function of a width of liquid layer h. The solution obtained by Sekiya et al. 
(2003) was for a completely-molten particle, so it corresponds to h/vp = 1 as 
pointed by a filled circle. The solutions for Kato et al. (2006, dashed curve) 
were for arbitrary values of h/rp from to 1 and Wmax — > for /i/rp — >• 0. 
We linearly interpolate from v^iax = at h/rp = to I'max = 0.112pfinrp//x 
Sit h/vp — 1 (solid line). 
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